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Abstract 



A proof of strong cosmic censorship is presented for a class of solutions 
of the Einstein-Maxwell equations, those with polarized Gowdy symmetry. 
A key element of the argument is the observation that by means of a 
suitable choice of variables the central equations in this problem can be 
written in a form where they are identical to the central equations for 
general (i.e. non-polarized) vacuum Gowdy spacetimes. Using this it is 
seen that the deep results of Ringstrom on strong cosmic censorship in the 
vacuum case have implications for the Einstein-Maxwell case. Working 
out the geometrical meaning of these analytical results leads to the main 
conclusion. 

O ' 

On 

^ ! 1 Introduction 



The Gowdy spacetimes are a class of solutions of the vacuum Einstein equa- 
tions with two commuting spacelike Killing vectors. Factoring by the symmetry 
leads to effective equations in one time and one space dimension. To date the 
Gowdy spacetimes are one of the best model systems for obtaining a mathemat- 
ical understanding of the dynamics of inhomogeneous solutions of the Einstein 
equations. It is known that close analogues of the Gowdy spacetimes exist in 
the case that the vacuum Einstein equations are replaced by the source-free 
Einstein-Maxwell equations (see for instance 4J). It is thus natural to ask to 
what extent the known results on the physical properties of Gowdy spacetimes 
can be generalized to the case with an electromagnetic field. This question is 
the subject of this paper. When referring to Gowdy spacetimes in the following 
the spatial topology is always assumed to be T 3 (a three-dimensional torus) 
except where the contrary is stated explicitly. 

A variety of deep mathematical results have been obtained concerning Gowdy 
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spacetimes, culminating in the work of Hans Ringstrom [TB], [T7], (T^]. In par- 
ticular it was shown in 19J that strong cosmic censorship holds in the class of 
Gowdy spacetimes. This means that for generic initial data within this class the 
corresponding maximal globally hyperbolic development is inextendible. (Pre- 
cise definitions of these terms are given later in the paper.) A relatively simple 
special case of the Gowdy spacetimes are those which are polarized. In that case 
the main field equation becomes linear and this simplifies many things. In the 
case of the Einstein-Maxwell equations it is also possible to make a definition 
of polarized Gowdy symmetry but in that case the main field equations are no 
longer linear. In fact, as shown explicitly below, when expressed in suitable 
variables they are identical to the full (non-polarized) Gowdy equations of the 
vacuum case. Hence the results of Ringstrom concerning the solutions of the 
partial differential equations can be taken over directly. On the other hand, the 
geometrical and physical meaning of the variables is different in the two cases 
so that the interpretation of the results is different. This will be worked out in 
detail in the following. It leads in particular to a result on strong cosmic cen- 
sorship for solutions of the Einstein-Maxwell equations with polarized Gowdy 
symmetry (Theorem [2]) . 

To put this result into context it is worth recalling for which matter mod- 
els strong cosmic censorship has been proved for any inhomogeneous solutions 
of the Einstein-matter equations. In [7] strong cosmic censorship was proved 
for solutions of the Einstein- Vlasov system with T 2 symmetry. The genericity 
assumption made there excludes the vacuum case. The matter content is used 
in an essential way. Although it is plausible that the methods of [7] could be 
extended to many other matter models they do not suffice to treat the case 
of a Maxwell field alone. Similar results have been proved for solutions of the 
Einstein- Vlasov system with spherical symmetry and hyperbolic symmetry in 
[8] . The results of [7] were extended to the case with positive cosmological con- 
stant in [20] . Strong cosmic censorship has been proved for the plane symmetric 
Einstein-scalar field system in [5T] . For those Gowdy spacetimes (not considered 
further in the following) with topologies other than T 3 cosmic censorship has 
been studied in the polarized case in [5]. The results which have been listed 
are those existing positive results known to the authors on strong cosmic cen- 
sorship of this level of generality. They are valid for inhomogeneous solutions 
of the field equations, with or without matter. They concern solutions which 
belong to some symmetry class but do not require smallness assumptions on the 
initial data. 

Apart from cosmic censorship, more detailed information is obtained about 
the qualitative properties of the solutions of the Einstein-Maxwell equations near 
the initial singularity and at late times. For generic solutions it is shown that the 
Kretschmann scalar R a p 7 sR a P 1 ' s blows up uniformly on certain regions as the 
singularity is approached. It is also shown that the spacetimes are geodesically 
complete towards the future and that they exhibit oscillatory behaviour at late 
times analogous to that found in [16] for Gowdy spacetimes. 

In Section [2] definitions of Gowdy and polarized Gowdy symmetry are in- 
troduced. These definitions, which make use of discrete symmetries, have not 
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previously been worked out in the literature and so are treated in detail here. 
In Section [3l a number of equations are displayed. They are the basic equations 
of the class of spacetimes which are of interest in the sequel. In Sections 0] and \5\ 
the main results on the asymptotic behaviour of solutions in the past and future 
direction respectively are proved. These results are combined in Section [5] to 
give the main result (Theorem [5]) and some thoughts are presented on how the 
results of this paper might be generalized in the future. 

This paper is based in part on the diploma thesis of the first author [12j . 

2 Gowdy spacetimes and discrete symmetries 

Let g a p be a Lorentzian metric on a four-dimensional manifold M and suppose 
that g a p is a solution of the Einstein equations coupled to some matter model. 
A symmetry of this solution is an isometry of the metric which also preserves 
the matter fields. The solution is said to have T 2 (or U(l) x U(l)) symmetry 
if this group acts on M by symmetries of the solution. Think of a solution 
as being defined by Cauchy data. If the Cauchy data are invariant under a 
symmetry group and the Einstein equations coupled to the given matter model 
have a well-posed initial value problem then there is an associated action of the 
group on the solution (cf. the discussion in section 9.2 of [T5]). 

Now suppose a spacetime is given where the metric takes the form 

Here 9, x and y are periodic coordinates on T 3 and the functions rj, U, A, G, H 
and R depend only on t and 9. This metric has T 2 symmetry with the group 
action generated by the Killing vectors ^ and J^. This is a form of the metric 
used in [2] . In that paper it was shown that in any globally hyperbolic spacetime 
with T 2 symmetry and a Cauchy hypersurface diffeomorphic to T 3 coordinates 
of this type can be introduced in a neighbourhood of any Cauchy hypersurface. 

A metric of the form |T]) will be said to have Gowdy symmetry if the trans- 
formation which simultaneously maps x to — x and y to —y is an isometry. It 
follows that G and H must vanish. Conversely, if G and H vanish the space- 
time has this extra symmetry. The terminology comes from the fact that the 
solutions of the vacuum Einstein equations with this symmetry are precisely the 
Gowdy spacetimes. A metric with Gowdy symmetry is said to be polarized if 
the individual transformations mapping x to —x and y to — y are symmetries. 
This is equivalent to the vanishing of A. Next a more abstract characteriza- 
tion in terms of discrete symmetries will be derived. Let G\ be the group of 
isometries of R 2 generated by translations and reflection in the origin. Let the 
latter transformation be denoted by R\^- G\ is a semidirect product of R 2 and 
Z 2 ■ Suppose that in a spacetime with T 2 symmetry the action of T 2 extends to 
an action of G\ by isometries. Then the coordinate form of the transformation 
corresponding to R\_2 is given up to an overall translation by the transformation 
which simultaneously sends x to —x and y to —y. This follows from the relation 
Ri,2Ta,b = T- a -bRi.2 which holds for any translation T a ,b- To see this, note 
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first that the isometry corresponding to R12 cannot be a translation. Hence it 
has a fixed point. Evaluating the relation at that fixed point gives the desired 
conclusion for the restriction to one orbit. Consider the geodesic in S through 
that point orthogonal to the orbit. Its intersection with any other orbit must 
also be a fixed point of R\2- Thus the necessary translation is the same for all 
orbits on the initial hypersurface. The same conclusion can then be obtained 
for any orbit in the spacetime by a similar argument using a geodesic of the 
spacetime through the fixed point which is orthogonal to the initial hypersur- 
face. In this way an abstract characterization of Gowdy symmetry in terms of 
an action of the group G\ is obtained. Similarly, polarized Gowdy symmetry 
can be characterized by the fact that the T 2 action extends to an action by 
isometries of the group G2 of isometries of R 2 generated by translations and 
the reflections Ri and R2 in the two Cartesian coordinates. The group G2 is 
a semidirect product of R 2 with Z2 ® Z%- The elements R\ and R2 satisfy 
the relations RiT a ,b — T- a ,bRi and i?2?a,fc = T a ,-bR2- It follows by the same 
type of argument as in the previous case that R± and R2 can be identified with 
reflections in x and y. 

Let us suppose for simplicity that the torus has unit coordinate volume with 
respect to the coordinates (x,y). Then the area of an orbit of the symmetry 
group contained in the domain of the above coordinate system is given by R. 
This can be used to extend the definition of R to the whole spacetime by re- 
quiring its value at a given point to be equal to the area of the orbit on which 
that point lies. The function R is called the area radius. It will be assumed 
in the following that the Maxwell field does not vanish identically so that the 
spacetime is not flat. Recall that the Maxwell field satisfies the dominant energy 
condition. A proof of this is given in the Appendix. It follows from Proposition 
3.1 of J3] that the gradient of R is always timelike. As a consequence the func- 
tion R can be used as a time coordinate. This is referred to as an areal time 
coordinate. It is not a priori clear whether in a given solution there is any level 
surface of R which is a Cauchy surface. An argument which shows that there 
is a level surface of this kind is presented in the next section. 

Having defined Gowdy and polarized Gowdy symmetry for metrics the next 
step is to do the same for solutions of the Einstein-Maxwell equations. To do 
this the action of the groups G\ and G2 must be extended to the Maxwell field. 
Any diffeomorphism transports any tensor field in a natural way. If the diffeo- 
morphism is an isometry consistency with the Einstein equations requires that 
the energy-momentum tensor must be transported in this way. It is, however, 
possible to transport matter fields in a different way and this will be done here. 
Consider the case that the electromagnetic field is defined by a four-potential 
A a . If <j) is a diffeomorphism consider the transformation (g, A) 1— > (4>*g, — </>*A). 
This leads to the transformation F 1— > —<fi*F for the electromagnetic field. This 
ensures that the energy- momentum tensor is left invariant. If the action of the 
group Gi is extended so as to transform the electromagnetic field in this way 
then Aq and A\ vanish. Denote the remaining components of the potential by 
u> = A2 and x = A3. This provides the definition of Gowdy symmetry for the 
Einstein-Maxwell system. The only non-trivial components of the electromag- 
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netic field are F 02 , F Q3 , Fi 2 and F13. Taking instead Gi provides the definition 
of polarized Gowdy symmetry. This does not lead to a direct further restriction 
on the potential but does further restrict the metric. It will be shown in the 
next section that the Einstein-Maxwell equations lead to a consistency condition 
which apparently forces one of the components of the potential to be set to zero 
in the polarized case. 

It would have been possible to use the more straightforward transformation 
(g,A) i— > (<j)*9, 4>*A) instead of that just introduced. This leads to a different 
class of solutions of the Einstein-Maxwell equations. In that case the only non- 
trivial components of the potential are Aq and A\. The only component of the 
electromagnetic field which may be non- vanishing in this case is then Fqi- If 
fields are allowed which do not arise from a potential then F23 may also be 
non-vanishing. Solutions of the Einstein-Maxwell equations of this type were 
considered in [22] . It seems that they exhibit oscillatory behaviour near the 
singularity and are thus more complicated than the solutions considered in the 
following. This other symmetry class is not considered further in the present 
work. 

3 Basic equations 

Consider Cauchy data for a solution of the Einstein-Maxwell equations with 
polarized Gowdy symmetry and suppose that these data have a constant value 
of R. In the solution evolving from these data coordinates exist on a neigh- 
bourhood of the initial hypersurface where the metric takes the form |T]) with 
G = H = A = 0. The function R can be introduced as a new time coordinate 
and after that a new spatial coordinate can be introduced so that goi — 0. Then 
the transformed metric is of the form (P} with G = H = A = and R = t except 
for the fact that dt 2 may have been replaced by an expression of the form adt 2 
for some function a(t, 0). Using the fact that for the Maxwell field T$ = — T/ 
one of the field equations (equation (11.42) on p. 200 of [TS]) implies that a t = 0. 
It is possible, while maintaining the coordinate conditions imposed up to now, 
to reparametrize the coordinate 9 on the initial hypersurface. Because a does 
not depend on time this can be used to impose the condition a = 1. Introduce 
new variables as follows: A = 4(?y — U) + logt, P = 2U — logt. Then the metric 
takes the form 

t~h?{-dt 2 + dd 2 ) + t(e p dx 2 + e- p dy 2 ) (2) 

In the new coordinates the original Cauchy hypersurface is not a hypersurface 
of constant t. That there does exist a Cauchy hypersurface of constant t can 
be shown as follows. The maximal Cauchy development of the initial data 
corresponds to a certain subset of (t, 6>)-space. Suppose that this set has a past 
boundary point p at which t is strictly positive. By using standard light cone 
estimates for the Gowdy equations it can be shown that P, x an d all their 
derivatives all orders are uniformly bounded on the intersection of the region 
corresponding to the maximal Cauchy development with a small neighbourhood 
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of p. By uniform continuity they extend continuously to p. Using arguments as 
in sections 4 and 5 of pQ or as in [3] it can be concluded that the solution can 
be extended to a larger region in such a way as to define a globally hyperbolic 
extension. This contradicts the definition of the maximal Cauchy development. 
Hence no such boundary point p exists. The solution exists all the way to t = 
and hence if to is any positive number smaller than the minimum of t on the 
original Cauchy hypersurface then t = to will be a Cauchy surface of the type 
whose existence was to be demonstrated. 

As explained in the previous section the only non-vanishing components of 
the four-potential defining the Maxwell field are — ui and A3 — \. With 
these choices the Einstein-Maxwell equations can be computed. The details are 
given in [12j . 



The full set of Einstein equations includes one involving the second derivatives 
of A. It is a consequence of (j6]) and ([7|) and plays a minor role in what follows. 
It is written down in terms of other variables later. The equation ([5]) is the 
consistency condition alluded to in the last section. The only large classes of 
solutions known which are consistent with this equation are those where either 
u> or x vanishes. In [4] the authors considered other possibilities but this did 
not lead to the discovery of other interesting classes of solutions. In most of 
what follows it will be assumed that u = 0. Assuming \ — would lead to an 
equivalent problem. These two cases are related by interchanging the roles of 
the two Killing vectors. 

The equations ©, ©, © and ® with u = can be shown to be equivalent 
to the equations for (not necessarily polarized) Gowdy vacuum spacetimes by 
introducing new variables. These are given by 



Then the equations for P, x and A coming from ©, (O, © and ([7]) are the 
basic equations for Gowdy vacuum spacetimes (cf. ([H]) with P replaced by P, 
Q replaced by x an d ^ replaced by A. These equations are: 



-P tt - t~ l P t + Pee = V P (-^ 2 + w |) - e p (- x 2 t + xl)] (3) 

~u> tt + wee = -Pt^t + Pe^e (4) 

-Xtt + Xee = - [-PtXt + PeXe] (5) 

At - t(P? + Pf) + 4[ e - p K 2 + u? d ) + e p ( X 2 t + xl)] (6) 

\g = 2tP t P 9 + 8(e- p cj t cj e + e p X tXe) (7) 

^tXt = ueX9 (8) 



P=-(P- log*) 
A=i(A-logi)-P 



(10) 



0) 



-P a - t~ l Pt + Pee = e 2p (~Xt + xl) 
-Xtt - t^xt + Xee = -2(-P t xt + PeXe) 



(11) 

(12) 
(13) 



G 



\g = 2t[P t P g + e 2p X tX0] 



(14) 



In terms of these variables the metric is 




(15) 



(16) 



4 Structure of the singularity 

Since the equations for solutions of the Einstein-Maxwell system with polar- 
ized Gowdy symmetry are identical to the well-known Gowdy equations in the 
vacuum case it can be concluded from a theorem of [10] that the solution cor- 
responding to any initial data with constant R exists globally on the interval 
(0, oo) when expressed in areal time. The purpose of this section is to investigate 
the asymptotics of these solutions in the limit t — ► 0. Interpreting a result of [T7] 
in terms of the present problem shows that for a fixed value of 9 the quantity 
t[(P 2 + e 2P xi )(*> 0)] 1 converges to a limit Vao(9) as t — > 0. This quantity is 
called the asymptotic velocity at 9. It turns out to be useful to examine regions 
Vg defined as the set of points (t,9) such that \9 — 9 Q \ < t. Intuitively this 
can be thought of as the causal future of a point of the singularity defined by 



The structure of the singularity will be described in detail for a generic set 
Q c of initial data, which is strictly analogous to a similar set introduced in [17] . 
using the relation between vacuum solutions and polarized Einstein-Maxwell 
solutions. By analogy with [T7| we can define what is meant by a solution 
having a non-degenerate false spike or a non-degenerate true spike at some 
value of 9. The detailed definition will not be reproduced here. For convenience 
call a point 9 an ordinary point for a given solution if < v ao (9) < 1 and 
lim t _»o tP t (t, 9) — —v 00 (9). Define Qi^ rn ,c to be the set of initial data determining 
solutions for which all points 9 arc ordinary except for I non-degenerate false 
spikes and m non-degenerate true spikes. Then Q c is the union of Gi, m ,c over all 
non-negative integers I and m. Data belonging to Q c are called generic. This 
subset is open in the set of initial data in the C° x C 1 topology and dense in 
the C°° topology [19]. A comment should be made on the subscript c which 
is included here for consistency with Ringstrom's notation. It is a reminder 
that initial data for the Gowdy equations must satisfy an integral constraint for 
topological reasons. An analogous constraint occurs in the Einstein-Maxwell 
case. It arises by integrating equation (TT4"]) with respect to 9. The integral of 
the left hand side vanishes due to periodicity and so the integral of the right 
hand side must vanish too. If there is a false spike at 9 then < v OCJ (9) < 1 while 
if there is a true spike 1 < Voo(9) < 2. Thus, in particular, generic data are such 
that Voc is never equal to one at any point. In the vacuum case this statement 
is enough to imply blow-up of the Kretschmann scalar along any past-directed 
causal curve due to Proposition 1.19 of [T7] . 



9 



7 



To prove inextcndibility it will be necessary to consider the case of normal 
points, false spikes and true spikes separately. First the case of a normal point 
will be examined. The expression for the Kretschmann scalar is unwieldy and it 
might be thought that it would be easier to prove inextcndibility in the electro- 
magnetic case by showing that the invariant F a pF al} blows up. Unfortunately 
this turns out not to be the case. There is a large class of Gowdy solutions 
(the low velocity solutions) which, when translated to the case of interest here, 
admit asymptotic expansions of the following form 

P(t,0) = -«oo(0) log t + ... (17) 

x (t,9) = q(e) + i/j(e)t 2v ~ + ■■■ (is) 

where the terms omitted are lower order as t — > 0. It follows that e~ 2A is 
proportional to t~ 2v °° in leading order while 

xl-X 2 t =4-±vl c &™-^ + ... (19) 

Since 

F^F afj = 2e- 2 \ X 2 e ~ X?) (20) 

it follows that the information contained in the above expansions and the fact 
that < < 1 do not always suffice to determine whether the invariant blows 
up or not. For if Voo — ^ and q 2 = Av^ij) 2 the leading term in the expansion 
for F a ^F a (i vanishes and whether the next to leading term blows up cannot 
be decided using the available information. For the solutions considered here 
F a P*F a p = and so this invariant is not helpful. 

Next the Kretschmann scalar will be examined. An observation which sim- 
plifies the calculations is that, due to the discrete symmetries of the spacetimes 
under consideration, a component R a f3-yS which does not vanish must have an 
even number of indices equal to two and an even number equal to three. It 
follows that up to symmetries (not including the algebraic Bianchi identity) the 
only components which may be non-zero are 

P0IOI7 -^0202, -^0212, -Rl212, -^0303, -^0313, -^1313, -^2323 (21) 

The Kretschmann scalar is equal to 

a ! p pOlOl , p p0202 , n pl212 , p p0303 
4[-Koioi-K + -K0202-K +-n-i212-n- + -K0303-K 

,p pl313 , p p2323l , orp p0212 1 p p0313i f Q\ 
+-n-1313-ft + -fl2323^l J + °[-K0212-fl + -K()313^ J- 

In terms of the metric coefficients the individual terms are given by 

t -4 e -4(A+P) M 2 Xtt _ t 2p tt + t 2 ( p g() + - Xg(j)] 2 
f -4 e -4(A + P }[ ^ t + t 2p~ Xt _ t 2p u _ t p t + t 2 { p2 + p^ 
t -4 e -4(X+P) [tXt + t 2~ Xt p t + t p t + t 2p2 _ t 2 { p ee _ p e - Xe)] 2 
t -4 e -4(X+P) [t 2p u _ t 2p~ Xt _ 2t 2p2 t 2 { p eXe + p2 )} 2 



0101 



R p0202 



-R0101-R 

?0202 

p pl212 
-K1212-K 

R R0303 
-KQ303-K 
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i- 4 e- 4 ^[-i 2 F t A t - t 2 P 2 + t 2 (P ee ~ 2P 2 - P e X e )] 2 

t -4 e -4(A+P) { _ t p t _ t 2p2 + t 2p2 ]2 

_ t -4 e -4(A+P) [ _^ _ f2 p t p g + t 2p te _ t 2p elt _ t 2 le p t] 2 
_ i -4 e -4(A + P) [ _ i 2 j p te + t 2p- Xe + t 2p glt + M 2p t p g] 2 

The first six terms are manifestly non-negative while the last two are manifestly 
non-positive. The fact that all terms contain a common factor helps when com- 
paring their magnitudes. It suffices to compare the magnitudes of the quantities 
in square brackets. Suppose it can be shown that in a given solution at least 
one of the positive terms tends to infinity as t — > and that both the negative 
terms are neglible with respect to that positive term. Then this implies that the 
Kretschmann scalar tends to infinity. It turns out, as will be explained in more 
detail later, that for this question the spatial derivatives in these expressions 
are unimportant compared to the time derivatives. The essential information 
needed for the proof is contained in the following expressions: 

t -A e -ACX + P) [ _ t 2- Xtt _ t 2p tt] 2 + 
t -4 e - 4 (A + P) [t ^ + t 2p~ Xt _ t 2p u _ t p t] 2 + 
t _4 e _4(X+P) [tSi + t 2 Xt p t + f p t + t 2p2 ]2 + _ 
t-4 e -4(A+P) [t 2p ti _ t 2p tXt _ 2t 2p2f + _ _ _ 

t- A e-< x+p )[-t 2 P t \ t -t 2 P?] 2 + ... 
t- i e-< x+p \-tP t -t 2 P?] 2 + ... 



The terms not explicitly written in these equations are not relevant for the proof 
which follows. It is not claimed that in the approach to the singularity each of 
the terms which have been omitted is small in comparison to the terms which 
have been written out explicitly. 

Consider now the case of a normal point Qq. Then Proposition 1.5 of [T7] 
provides information about the limiting behaviour of the solution in a neigh- 
bourhood of 6*o- It implies that tPt converges to —v^ and that tPg and t 2 Pte 
tend to zero. It also shows that if equations (fT3|) and (fT4|) are multiplied by t 
the terms containing derivatives of x on the right hand sides tend to zero as 
t — > 0. If follows that tX t converges to v 2 ^ and tXg tends to zero as t — > 0. These 
facts together show that the quantities in square brackets in the expressions 
for i?02i 2 -R 0212 and i?03i3-R 0313 tend to zero at the singularity. The expression 
in the square bracket in -R2323-R 2323 converges to Voo — v 2 ^ ^ 0. The expres- 
sion e _4 ' p+A - ) is asymptotic to an expression proportional to t 4 ( v °°~ v <*>h Since 
^oo(l— Woo) < 1 this is enough to conclude that the curvature blows up uniformly 
in a neighbourhood of a normal point 9q. It is also part of the conclusions of 
Proposition 1.5 of [17] that if 9q is a false spike a change of variables can be 
carried out which gives a new solution with the asymptotic properties which 



p pl313 
-n.131.3-n. 



tl313 J 
^2323^ 
^0212-i 
?0313 J 



p p2323 
-n-2323-tl 



R f?0212 
-K0212-K 



p p0313 



p pOlOl 

p p0202 
-TC0202-K 

1212 



R1212R 
^0303 
^1313 
^2323 
?0212 
?0313 



p p0303 
-nV13n3.fi 



p pl313 
-fll313-fl 



p p2323 
-fl 2 323-<l 



p p0212 
-nVj212-<l 



p p0313 
-nV)313-fl 
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hold for a normal point. Since it is known that this transformation does not 
change the spacetime geometry but only the way it is represented the curvature 
properties are unaffected. Thus curvature blow-up also holds in the case of a 
false spike. It remains to treat the case of a true spike at a point 0q. In that case 
Proposition 6.7 of [T7] implies that tPg and te p \e converge to zero uniformly 
on T>e . Corollary 6.8 shows that t 2 P t e also converges to zero there. By Propo- 
sition 6.11 of [17] and the remark following it tP t and te P Xt converge to — Woo 
and zero respectively on that region. In this case Woo(l — «oo) < 0. Thus the 
Kretschmann scalar can be controlled as in the other cases, with the difference 
that the control is only obtained on the region T>g rather than for an interval 
of 6 containing 9 . 

In the following theorem the topology on the set of initial data used is the 
C°° topology. 

Theorem 1 There is an open dense subset Q c of the set of smooth initial data 
for the Einstein- Maxwell equations with polarized Gowdy symmetry and con- 
stant areal time such that the Kretschmann scalar tends to infinity along any 
inextendible past-directed causal geodesic. 

Proof The set Q c has been defined above. For a solution evolving from data 
belonging to this set each point 0q is either an ordinary point, a false spike or a 
true spike. It has been shown that in each of these three cases the Kretschmann 
scalar blows up uniformly on T>g as t — > 0. It is shown in [17j that any in- 
extendible past-directed causal geodesic eventually lies in a region of this type 
with some value of 9o as t — > 0. The proof of this is pure Lorentzian geome- 
try and makes no use of the Einstein equations. Hence it also applies to the 
situation of this theorem. 

5 Future geodesic completeness 

This section is concerned with the asymptotic behaviour of solutions in the fu- 
ture, i.e. with the behaviour as the areal time t tends to infinity. The asymptotic 
behaviour of solutions of the vacuum Gowdy equations has been determined in 
detail in [16j . This information suffices in particular to prove future geodesic 
completeness in the vacuum case. It will be shown that similar arguments 
can be used to give a proof of future geodesic completeness of solutions of the 
Einstein-Maxwell equations with polarized Gowdy symmetry. 

As in the vacuum case it is convenient in the proof of geodesic complete- 
ness to distinguish between homogeneous and inhomogeneous spacetimes. The 
homogeneous case is that in which P, x and A depend only on t. Geodesic com- 
pleteness in the homogeneous case follows from Theorem 2.1 of [13]. That the- 
orem applies to solutions of the Einstein equations coupled to a general class of 
matter models which satisfy two conditions, the dominant energy condition and 
the non-negative sum pressures condition. Both of these hold for the Maxwell 
field. In fact the second follows from the first because the energy-momentum 
tensor of the Maxwell field has vanishing trace. It should also be noted that for 
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these solutions, which are of Bianchi type I, the mean curvature of the hyper- 
surfaces of constant t is everywhere negative and hence, in particular, bounded 
above. 

It remains to consider the inhomogeneous case. From Proposition 1.8 of 
[TB] it can be concluded that P t , Pe, e p \t and e P xe are 0{t~ x / 2 ) as t — > oo, 
uniformly in 9. This means in particular that At is bounded in the future. From 
Theorem 1.7 of [16] it follows that A = c^i + (9(log t) for a positive constant c\ as 
t — > oo. It is at this point that the condition that the solution is inhomogeneous 
is used. Combining these facts shows that ' < A + P = c^t + Oft?). Define an 
orthonormal frame {e^} by normalizing the coordinate basis vectors Jj, g|, 
j^. A computation shows that the corresponding rotation coefficients satisfy 
an inequality of the form 

\g{V eit e v ,e p )\<Ct-%e- x - p (23) 

for a constant C whenever at least two of the indices (/x, v 1 p) are equal to two 
or three. The facts just listed suffice to show that the argument used to prove 
geodesic completeness in the vacuum case in [IB] also applies to the Einstein- 
Maxwell case, proving geodesic completeness in that case too. 

In [16] further information about the late-time behaviour of Gowdy space- 
times is derived. It is shown that, remarkably, there are solutions with os- 
cillations of P which persist for all times. More precisely, the amplitude of 
the oscillation does not decay and the fact that P t converges to zero is ex- 
plained by the fact that the frequency of the oscillation decreases with time. 
The derivative Pe does decay, which means that in a sense the solution becomes 
homogeneous asymptotically. Nevertheless the dynamics of these solutions is 
different from that of any homogeneous solution. These solutions can be used 
to produce solutions of the Einstein-Maxwell equations by the correspondence 
which is central to this paper. The statements made about P in the vacuum 
case apply directly to P in the electromagnetic case. Considering \ it is seen 
that the electromagnetic potential oscillates with an amplitude which does not 
decay while the decrease in the frequency of the oscillation leads to decay of the 
electric field. The potential becomes spatially homogeneous asymptotically and 
this gives decay of the magnetic field. 

A key step in proving the results about future asymptotics in the vacuum 
case is to show that the following energy- like quantity tends to zero as t — > oo. 

£(*) = \ I \Pt + Pe + z 2P (Q 2 t + Qe)](t, 6)d0. (24) 

Of course there is a corresponding quantity in the Einstein-Maxwell case which 
also tends to zero as t — > oo. It is given explicitly by 

Ht) - ~ f [Pt + Pe + e 2P (x 2 t + Xlm 0)dS. (25) 
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6 Conclusions and outlook 



Given a matter model in general relativity for which the Einstein-matter equa- 
tions have a well-posed initial value problem there exists a maximal Cauchy 
development corresponding to each initial data set. The strong cosmic censor- 
ship hypothesis says that for generic initial data this maximal development is 
not extendible to a larger spacetime. There are some technical choices involved 
in this definition: how much regularity is assumed for the data, what topology 
is put on these data and what regularity properties are assumed of an extension. 
Here no general discussion of these issues will be given - we will just make one 
reasonable choice. It is assumed that the data are smooth, i.e. C°°. They are 
topologized using the C°° topology, i.e. the topology of uniform convergence 
of a geometrical object together with its derivatives of all orders. Extendibil- 
ity is taken in the sense of extensions of the geometry of class C 2 . It should 
be noted that since cosmic censorship is such a hard problem it is common to 
consider the related problem where the initial data is generic subject to some 
extra restrictions, e.g. symmetry restrictions. This is the case in the results of 
this paper. 

The main theorem of the paper can now be given. It is a statement about 
strong cosmic censorship for the class of spacetimes involved. 

Theorem 2 For data belonging to the open dense subset Q c of the set of smooth 
initial data for the Einstein- Maxwell equations with polarized Gowdy symmetry 
and constant area radius the corresponding maximal Cauchy development is in- 
extendible. 

Proof Let (g,A) be initial data belonging to Q c . By the results of Section 
0] the maximal Cauchy development of these data is inextendible towards the 
past since the Kretschmann scalar is unbounded along any inextendible past- 
directed causal geodesic. By the results of Section [S] the maximal development 
is future geodesically complete and hence inextendible towards the future. This 
completes the proof. 

In fact inextendibility towards the future already follows from global exis- 
tence to the future in areal time, as has been shown in [6]. It is nevertheless 
good that more has been proved about the future asymptotics since this gives 
more insight into what is going on. It should also be noted that the argument 
using global existence in areal time is fundamentally restricted to spacetimes 
with symmetry while proving geodesic completeness is a strategy which can in 
principle be applied to situations without any symmetry. 

Note that the fact that the above theorem deals only with initial data with 
constant area radius is not a serious restriction since any solution admits a 
Cauchy hypersurface carrying data of this type. 

In this paper the dynamics of solutions of the Einstein-Maxwell equations 
with polarized Gowdy symmetry has been analysed rather comprehensively. 
This was based on the existing analysis of the dynamics of vacuum Gowdy space- 
times which need not be polarized. This raises the obvious question, whether it 
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is possible to do something similar for general solutions of the Einstein-Maxwell 
equations with Gowdy symmetry. Some aspects of this problem will now be 
discussed. In the full Einstein-Maxwell case there are two metric coefficients P 
and Q as in the vacuum case and two components u> and \ of the potential. 

In looking at the behaviour of the solutions in the past a possible starting 
point is to try to extend the results on the existence of low velocity solutions of 
the Gowdy equations proved in [9] using Fuchsian methods. A more complicated 
system of this type, the Einstein-Maxwell-dilaton-axion system was discussed 
from this point of view in [TT] . Setting the dilaton, the axion and the dilaton 
coupling constant in that system to zero gives rise to the Einstein-Maxwell 
system. A special case considered in [TT] is in principle relevant to this but 
unfortunately it is not self-consistent in the non-polarized case since one of the 
field equations was neglected (cf. the remarks on p. 166 of [15]). Other special 
cases considered in [11] may provide an avenue to a greater understanding of 
singularities of solutions of the full Einstein-Maxwell-dilaton-axion system. In 
any case, the question of whether Fuchsian techniques can be used to obtain 
solutions of the Einstein-Maxwell equations with Gowdy symmetry depending 
on the full number of free functions (i.e. eight) remains open. 

Global existence in the future in areal time is known and this suffices to 
prove the part of strong cosmic censorship which concerns the future evolution 
using the arguments of [3] . On the other hand the detailed asymptotics have not 
been determined. Parts of the argument which applies to the vacuum Gowdy 
case have been extended to the Einstein-Maxwell case in [18]. These are the 
parts which concern the decay of an energy-like quantity analogous to those 
defined in ([24]) and (f25]l . However other important parts of what was done in 
the vacuum case have not yet been generalized. 

These remarks show that the task of extending the results on the Einstein- 
Maxwell equations with polarized Gowdy symmetry obtained in this paper to 
the case of solutions with fully general Gowdy symmetry presents some inter- 
esting challenges. 



Appendix The dominant energy condition for the Maxwell field. 

It is common knowledge that the Maxwell field satisfies the dominant energy 
condition but it is difficult to find a concise, elementary and self-contained proof 
in the literature. For the convenience of the reader one will be included here. 
The dominant energy condition for an energy-momentum tensor T a ^ can be 
formulated as the condition that T a pv a w^ > for all future-pointing causal 
vectors v a and w 13 . Since any future-pointing causal vector can be written as a 
sum of two future-pointing null vectors it suffices to verify the condition in the 
case that v a and w a are equal to future-pointing null vectors l a and n a . By 
possibly rescaling l a and adding two spacelike vectors x a and y a these may be 
completed to a null basis where the vectors have the following inner products: 
l a n a = -1, l a x a = l a y a = n a x a = n a y a = 0, x a x a = y a y a = 1, x a y a = 0. 
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The electromagnetic field can be expanded in this basis to give 

T aP = Al [a n f3] + Bl [a x p] + Cl[ a yp] + Dx [a y^ + En^x^ + Gn [a y p] (26) 
where the coefficients are constants. Straightforward computations give 

F ai F p ^rn p = J(A 2 + EB + GC), F a pF a/3 = ^(D 2 - A 2 ) -EB-GC (27) 

and 

W 3 = ^-(A 2 + D 2 ). (28) 
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